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Abstract 

Many types of point singularity have a topological index, or 'charge', associated with 
them. For example the phase of a complex field depending on two variables can either 
increase or decrease on making a clockwise circuit around a simple zero, enabling the zeros 
to be assigned charges of ±1. In random fields we can define a correlation function for the 
charge- weighted density of singularities. In many types of random fields, this correlation 
function satisfies an identity which shows that the singularities 'screen' each other perfectly: 
a positive singularity is surrounded by an excess of concentration of negatives which exactly 
cancel its charge, and vice- versa. This paper gives a simple and widely applicable derivation 
of this result. A counterexample where screening is incomplete is also exhibited. 
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1. Introduction 

This paper provides a simple and general explanation for a feature which has been 
noted in point singularities of several types of randomly defined functions. The phe- 
nomenon is most easily described in terms of a specific example. Consider the set of zeros 
of a complex valued random function <fi, depending on two real variables x = (x,y). It is 
assumed that the statistical properties of 0(x) are translationally invariant. The phase of <fi 
may either increase or decrease by 2tt on traversing a clockwise circuit about a simple zero. 
Accordingly, the zeros may be described as carrying either positive or negative charges. It 
has been noted that positive zeros tend to be surrounded by negative ones, and vice- versa: 
by analogy with models of ionic fluids and plasmas this has been described as a 'screening' 
effect. The effect is expressed quantitatively in terms of a correlation function C describing 
the charge- weighted density of zeros (a precise definition will be given in section 2). Perfect 
screening is characterised by the relation 



where the integral extends over the entire two-dimensional plane. 

A screening relation of this type was first noted by Halperin [1] as a consequence 
of an analytical evaluation of such a correlation function in the case of zeros of complex 
functions. The screening relation was also discussed by Liu and Mazenko [2] in a related 
context. Later, the effect was surmised to exist in degeneracies of a random matrix model 
used to investigate the quantised Hall effect [3]. Numerical demonstrations of perfect 
screening for both the zeros and extrema of Gaussian random fields in two dimensions 
have been published [4]. More recently, analytical evaluation of correlation functions for 
several types of point singularity has provided further evidence that (1.1) is valid in many 
cases [5-7]. No general explanation appears to have been published before now. 

The explanation given here is very simple and widely applicable, and leads one to 
expect that the perfect screening described by (1.1) may be present wherever there is a 
charge-neutral gas of point singularities. It depends upon an assumption about another 
correlation function, which is verified in a very general context for random fields with a 
Gaussian distribution, and which is expected to be much easier to establish than (1.1) 
in other cases. One exception to (1.1), relating to components of eigenvectors of random 
matrices, is described and explained. This counterexample is of some physical interest 
because it arises in the topological description of the integer quantised Hall effect, discussed 
in references [8-10]. 

This paper will use (A) to denote the ensemble average of a quantity A. Because of 
the assumption that the statistics are translationally invariant, the correlation function 
between 0(x) and 0(x') is a function of X = x — x' only. The correlation function may 
therefore be written (</>(x)0(x')) = c(x — x'). It is distinct from that occuring in (1.1), 
which describes the singularities of the field. 

2. A derivation of the screening relation 

Equation (1.1) will be derived for the case mentioned in the introduction, namely zeros 
of a complex function of two real variables, (f>(x, y), but the approach is easily generalised. 




(1.1) 
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The total charge enclosed by a circuit C is 
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ds F(s) (2.1) 



where ds is an element of distance around the boundary. Here F(s) is shorthand for the 
integrand: for other types of singularity where their number can be determined from a 
line integral F(s) would be replaced by a different function. It will be assumed that in the 
more general case F(s) is a 'local' function of 0, depending only on the function and its 
derivatives at s. The integral can only be different from zero if the region A enclosed by C 
contains at least one zero of 4>. The zeros are located at positions Xj. The total charge is 

Q= ft (2-2) 

{i|xie^l} 

where qi are the charges of the individual zeros, which take the values ±1. A charge- 
weighted density of zeros, p(x), is defined: 

p(x) = J]g i «J(x-x i ) . (2.3) 

i 

The total charge Q enclosed by C can also be expressed in terms of this density: 

Q = [ dA p(x) . (2.4) 
J A 

Because it was assumed that the random process generating </>(x) is translationally in- 
variant, correlation functions depend only upon the difference between coordinates. We 
consider random processes which are symmetric between positive and negative charges, so 
that (p(x)) = 0. We define the mean density of zeros (without charge weighting) to be po- 
The correlation function that we consider, already referred to in (1.1), is defined by 

C(X) = (p(x + X)p(x)> . (2.5) 

This correlation function is the sum of a singular part, po5(X), and a regular part, which 
according to (1.1) is expected to cancel the weight of the delta function when integrated 
over all X. 

Now consider the correlation function of the quantity F(s) introduced in (2.1). The 
correlation between F(s) and F(s') depends upon both the distance between the points 
labelled by s and s', and on the angle between the tangent vectors to the curve C at these 
points. Within broad limits the choice of the curve C is immaterial to the argument, so 
that this may be taken to be a circle of radius R. In this case: 

(F(s)F(s')) = D(s - s', R) = D (s - s') + 0(1/R) (2.6) 



where Dq is the limiting form of the correlation function which applies when the tangent 
vectors at s and s' are parallel. 

In order to prove (1.1) it will be assumed that the correlation functions Dq(s — s') 
decays faster than l/|s — s'\ as |s — s'| — > oo. If (j) is a Gaussian random function, it 
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is possible to show that this condition is satisfied if the magnitude of the correlation 
c(X) = (0(X + x)</>(x)) decreases sufficiently rapidly as |X| — > oo. The argument is 
described in an Appendix; it is applicable to other types of point singularities as well as 
zeros of complex functions. 

Using (2.1) and (2.6), (Q 2 ) is expressed in terms of the correlation function of F: 

(Q 2 ) = [ ds [ ds' D(s - s', R) . (2.7) 
Jc Jc 

Assuming that the correlation function Dq decays faster than l/\s — s'\ as \s — s'\ — > oo, 
in the limit where R is large, the integral (2.7) is dominated by the region where s — s' is 
small: 

/oo 
ds D (s) + 0(1) . (2.8) 
-oo 

To gain information about the correlation function C(X), consider the mean-squared 
charge: from (2.4) we have 

(Q 2 ) = [ dx I dx' C(x - x') . (2.9) 
J A J A 

Again, consider the case where A is the interior of a circle radius R. If the integral on the 
left hand side of (1.1) exists, then the integral in (2.9) is dominated by contributions from 
where x is close to x', so that 

(Q 2 ) = irR 2 J dX C(X) + O(R) . (2.10) 

Equations (2.8) and (2.10) are only consistent if the term proportional to R 2 in (2.8) 
vanishes. This implies that (1.1) is satisfied. This argument has established that (1.1) is 
correct for the case of zeros of a complex function. 

If <p is not Gaussian, there does not appear to be any general argument which indicates 
that D (s — s') decays sufficiently rapidly. The argument given in the appendix can be 
adapted to various types of no n- Gaussian field, but different examples must be treated 
on a case-by-case basis. However, it is expected to be much easier to establish that the 
integral (2.8) converges than to establish (1.1), which is an exact identity satisfied by a 
highly singular function of the underlying field (p. 



3. Generalisations 

The argument presented in section 2 may be extended to many different random fields 
in different dimensions: there is no requirement that it should be possible to determine 
the correlation function C(X) exactly, and the argument can be extended to non-Gaussian 
fields. 

The argument for perfect screening might be applicable to any point singularities 
which carry indices which we term 'charges', when the gas of charges is, on average, neutral. 
The critical requirement is that the charge within a closed region can be determined from 
a surface integral involving a field F which was derived in some way from the original 
random field (p. It was assumed that this secondary field has a correlation function which 
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decays sufficiently rapidly at infinity, and which depends only upon the separation of the 
pair of coordinates. 

Most charged point singularities can be detected by using a surface integral. As a 
second example, consider stationary points of a real function of two real variables. These 
may be characterised by the Poincare index, which is defined by taking a clockwise circuit 
around the singularity. The index is +1 if the angle of the gradient vector rotates by 
+2-7T (i.e., in the same direction as the circuit), and —1 if the gradient vector rotates 
in the opposite direction to the circuit. Thus maxima and minima have index +1, and 
saddles have index — 1. Perfect screening has been also been demonstrated in this case 
by calculating the correlation function exactly [1,2], but it is instructive to see how the 
argument of section 2 is adapted in this simple case. 

If x(s) is the closed curve C, and v(s) = V0(x(s)) is the gradient vector at s, then 
the total Poincare index for the stationary points within C is 



where the second equality defines F(s) in this example. All of the arguments that were 
applied to the function F(s) defined in (2.1) are equally valid for that defined in (3.1). It 
follows that the correlation function of the density of extrema weighted by their Poincare 
index also exhibits perfect screening under quite general conditions. 

The argument also extends directly to cases in higher dimensions, for example it 
can be used to explain the example of perfect screening which was discussed in reference 
[3]. The reasoning for this case will be summarised briefly: the reader should refer to 
the earler papers for definitions of the quantities. Reference [3] considered the density of 
degeneracies between pairs of levels in random Hermitean matrices, which were a periodic 
function of three parameters, xi, x^ and X3. There is a topological charge, the Chern 
index, associated with the energy levels of a two-parameter family of Hermitean matrices 
which are periodic in two parameters, x\ and xi- The Chern index is an integer which may 
represent a quantised Hall conductance associated with each energy level [8,9]. Varying the 
third parameter xs allows pairs of levels to become degenerate. When two levels become 
degenerate, the Chern index on one level increases by one, while that of the other level 
decreases by one. This enables the degeneracies to be assigned a charge of ±1, depending 
upon the sign of the change of the Chern index of the upper level resulting from increasing 
Xs [10]. Reference [10] also shows that the total charge of the degeneracies of a given level 
(with label n, say) within a closed region of the three parameter space is equal to the 
integral of the Berry phase two- form V n [11] over the surface of the region. 

The random matrix model discussed in reference [3] has statistical properties which 
are translationally invariant in x = (#1, #2, #3) space, and this symmtetry means that 
degeneracies are equally likely to have either sign. This is analogous to the situation 
described in section 2: we have a homogeneous distribution of 'particles' (degeneracies) 
which are equally likely to have positive and negative charges, and the total charge Q within 
a three dimensional region is obtained by integrating a function, the Berry phase two-form, 
over its surface. Considereing a spherical volume of radius R, and writing {Q 2 ) in terms 
of the correlation function of charge density, gives an expression analogous to (2.8), with 
the leading term being 47ri? 3 /3 multiplied by the integral of the charge density correlation 
function. Expressing (Q 2 ) in terms of the surface integral leads to an expression analogous 
to (2.10) in which the leading term is AtxR 2 multiplied by an integral of the correlation 




(3.1) 
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function of the two-form, C(x — x') = (V r n (x)V r n (x / )). Provided this correlation function 
vanishes faster than l/|x — x'| 2 , that integral converges, and the correlation function of 
the charge density satisfies (1.1) (with the integral now being evaluated over a three- 
dimensional space). This argument explains the screening relation noted in reference [3]. 

4. A counterexample: zeros of eigenvector components 

The discussion in the previous section indicates that the argument explaining perfect 
screening is very general, and it might be suspected that perfect screening is universal in 
charge-neutral gases of singularities. However this section will describe a counterexample, 
which arises from the same random matrix model as was considered in reference [3]. The 
mathematical structure of the model will be explained, but the reader should refer to [3] 
for a discussion of the physical motivation of this model. 

Consider a complex Hermitean N xN matrix H, with elements H nm (xi, X2) which are 
a function of two real parameters, x\ and x^- The matrix elements H nm can be Gaussian 
random functions, depending smoothly on x = (x\, x%), with a correlation function between 
H nm (x) and H n i m i(-x!) which depends only upon x — x' and which decays faster than 
l/|x — x'| as |x — x'| — > 00. In this case each component of any eigenvector of H is a 
complex random function of x\ and X2- We consider one eigenvector component which will 
be termed (f>(xi,X2)- The function <f)(xi,X2) is regular at (xi,X2) unless the corresponding 
eigenvalue is degenerate at (xi, X2). Three parameters must be varied to cause degeneracies 
of eigenvalues of a complex Hermitean matrix, so that we expect <f>(xi,X2) to be regular 
everywhere in the (^1,^2) plane. 

Let us consider the zeros of <fi. There is nothing to favour one index of these zeros 
over the other, so that the distribution of zeros across the (xi,X2) plane is expected to 
be charge neutral. We can ask whether the screening relation applies to the zeros of the 
eigenvector component (f>(xi,X2). It will demonstrated that in the case of 2 x 2 matrices 
screeening the screening is not perfect. It will be argued that screening is also not perfect 
when N > 2. 

The indices of zeros of a component of an eigenvector of a two-parameter familiy of her- 
mitean matrices play a role in the topological characterisation of the integer quantised Hall 
effect [9] . In the case of independent electrons moving through a perfect two-dimensional 
crystal, the electron states form bands labelled by two Bloch wavevectors, x\ and x^.- The 
Hamiltonian is a periodic function of x\ and X2, with a unit cell which is termed the 
Brillouin zone. The Hall conductance of a band is Qe 2 /h, where Q is an integer. One 
way to calculate Q involves looking at any component of the eigenvector defining the 
wavefunction of the band. All of the zeros of within a Brilloiun zone are located, and 
thier indices qi = ±1 are determined. The Hall conductance integer Q is then the sum of 
the indices: Q — ^ t^. 

In cases where the electrons move in a simple periodic potential the integers Q can 
be large (although these situations would be very hard to probe experimentally). If the 
system is disordered, or the unit cell is large, it is reasonable to propose using a random 
matrix model for the statistical properties of the Chern numbers. An appropriate random 
matrix model is described in [3] . The mean Chern number must (by symmetry) be equal to 
zero: (Q) = 0. The simplest statistical characterisation of the Chern numbers is through 
their variance (Q 2 ). 

The random matrix model in [3] was investigated by a combination of analytical and 
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numerical approaches. Numerical experiments reported there support an expression of the 
form 

(Q 2 ) = aAa 2 p 2 (4.1) 



where 



det 



/ dE n dE n \ 
\ dxi dxj I 



(4.2) 



is a measure of the sensitivity of energy levels to perurbation, p is the density of states, 
A is the area of the Brillouin zone, and a is a constant, determined numerically to be 
approximately 0.2 in the limit where the dimension of the matrix is large. Because equation 
(4.1) is proportional to A, perfect screening does not apply in the case of zeros of eigenvector 
components. We have already seen that screening is predicted by an argument with a very 
broad range of applicability. It is desirable to understand why screening fails at two levels. 
Firstly, which of the criteria stated in the derivation are not met? And secondly, can it be 
seen by an explicit calculation that screening is absent? 

First, let us consider the reason why the demonstration presented in section 2 cannot 
be applied. The argument uses a function F(s) which yields the phase change upon 
integration around the boundary, and the function which was defined there, F(s) = 
(l/0)lm(d0/ds), is also appropriate in this problem. It was assumed that the function 
F(s) has a correlation function which is statistically stationary (such that on the circular 
boundary (F(s)F(s')} is a function of s — s' only, and that this correlation function decays 
faster than l/\s — s'\. In the case where <p is a Gauss random function with a specified 
correlation function which is statistically stationary and rapidly decreasing, these assump- 
tions are certainly valid. If is a component of an eigenvector of a random matrix, these 
assumptions about F(s) may be challenged. In order to apply (2.1) we must asume that 
(f)(x\,X2) is a smooth function of x\ and xi- The eigenvectors of a matrix may be multi- 
plied by any complex number of modulus unity, exp[i#(:ri, x^f\ and remain eigenvectors. 
When constructing <p(xi,X2), from the random matrix H(xi,X2), the eigenvectors can be 
produced by some fixed algorithm (for matrices of dimension N > 4, this is necessarily 
numerical). This results in a 'raw' form for <l>(xi, X2) which is both statistically stationary 
and locally correlated. However, this function will have discontinuities in phase which 
must be removed by multiplying by a factor exp[i9(xi, X2)] (the function 0(xi,X2) is not 
unique). The function 6(xi,X2) cannot be constructed by any locally defined algorithm. 
It is therefore possible that the resulting regularised <f>(xi,X2) may be non-stationary, or 
non-locally correlated, or both. 

This general argument is a little unsatisfying, because it does not show that the perfect 
screening effect must fail. However the source of the failure can be seen clearly in the case 
of a 2 x 2 hermitean random matrix, with elements 



where x = (x\,X2), and we assume that the real- valued random functions /j(x) satisfy 




(4.3) 



</i(x)> = 



(4.4) 



(/,(x + xo)/,(x )) =<%C,(|x|) . 



(4.5) 
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The eigenvalues E± and corresponding normalised eigenvectors u-t- are 



E± = ±y/fi + f$ + fi (4.6) 

u ± = , eXp[ ^ ±] ( £ + ^ ) (4.7) 

where all the variables are real functions of x, and where 6>±(x) is chosen so that the 
first component of u± is a regular function of x. Let us consider the first component 
of the eigenvector corresponding to the E + branch, calling this function 0(x). Writing 
V>(x) = / 2 (x) + i/ 3 (x), we have 

0(x) = 6XPW = . (4.8) 

Consider the behaviour of cj> at a zero of ip. At each zero of ip, we have E + = \fi\. We must 
consider two cases, depending upon whether f\ is positive or negative. If fi is negative at 
the zero of ip, then in the neighbourhood of this zero we have 

<j> (4.9) 

so that (f) has a zero with the same index as ip, 9 having no singularity at these points. In 
the case where fi is positive at the zero of ip, in the neighbourhood of this point we have 

exp[ifl]-0 

At these zeros of ip, <p does not have a zero, and # must have a singularity in which it 
increments by ±2tt on making a circuit around the zero of ip, in order to cancel the phase 
singularity of 

We can now give a clear picture of why screening is not perfect in the eigenvector 
component <p, for the special case of this 2x2 random matrix model. We have seen that 
the zeros of <fi have the same indices as those of if), which do exhibit perfect screening. 
However, only a randomly chosen half of the zeros of ip (selected by the criterion that fi 
is negative at the zero) are represented as zeros of (p. It is therefore not expected that the 
delicate balance implied by perfect screening will be present in the zeros of cf>. In the limit 
where the correlation length of /i(x) is made short compared to that of / 2 and fs, the 
deleted zeros are selected completely randomly. 

This conclusion can be expressed quantitatively as follows. The sign-weighted density 
of zeros of the eigenvector component 0(x) is 

/V(x)=x(/i(x))^(x) (4.11) 

where p^,(x) is the sign-weighted density of zeros of ip = ji + i/ 3 , and the factor x(fi( x )) 
selects those zeros for which /i(x) is negative: x( x ) * s unity if x is negative, zero otherwise. 
Using the fact that fi is independent of ip, we have 

C (x) = (x(A(x ))x(A(xo +x))>CV(x) 
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= C x (x)CV(x) (4.12) 

where the second equality defines C x (x). The function satisfies equation (1.1). The 
function need not. For example, if fi has a correlation length which is short compared 
to /2 and /3, the correlation length of C x will also be short compared to that of C^. 
Equation (4.12) shows that in this limiting case the screening effect would be absent from 
the eigenvector component <p. 
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Appendix 

This appendix discusses the calculation of the correlation function (F(s)F(s')), in the 
limit where the separation \s — s'\ is large. It is assumed that F(s) is a function of the field 
4>(s) and its derivative (j)'(s), that is F(s) = J-((j)(s),(j)'(s)) for some function !F(a,b), but 
the argument is easily extended to cases where F(s) depends on higher derivatives. The 
argument can also be adapted to higher dimensions. The calculation discussed here does 
use the assumption that the field 4> is Gaussian, but extensions to non-Gaussian fields are 
possible. 

It will be shown that (F(s)F(s')) decays no more slowly (as \s — s'\ — > oo) than 
the most slowly decaying of the correlations (0(s)0(s')), (0 / (s)0'(s / )), ((/)' (s)0(s')) . The 
correlation functions of derivatives are related to the derivatives of the correlation function, 
for example if (0(s)0(s')> = c(s - s'), then (0 / (s)0 / (s / )) = ~c"(s - s'). 

The joint probability density for N Gaussian random variables (xi, X2, ..,xn) t = X, 
all of which satisfy (xi) = 0, is 

A*(X) = [(27r) 7V det(7]" 1/2 exp(-iX T (7- 1 X) (A.l) 

where C is the correlation matrix of the random variables x^, with elements 

Cij = (xiXj) . (A.2) 

Now write X = (x!,x 2 ) T , xi = (0(s), (p'(s)) T , x 2 = (<j>(s'), 4>'{s')) T , and write the correla- 
tion matrix (A.2) in the form 

r_( A e(s - s')a(s - s')\ . , 

U ~ \e(s - s')a(s - s') A ) 1 6) 

where A and a are 2x2 matrices; A is independent of s — s' , a has elements which are 
bounded as \s — s'\ — > oo, and e is a function of s — s' which decays no more rapidly than 
the most slowly decreasing correlation function of <p and <f>' as \s — s'\ — > oo. 
When |e| is sufficiently small, the inverse of C is approximated by 

C 1 = [_ eA -i~ aA -i A -i J+0(e 2 ). (AA) 

We also have det(C) = [det(i)] 2 + 0(e 2 ). Using (A.4) and (A.l), the joint probability 
density for xi and x 2 is 

^( Xl ' X2 ^ = Tn 7Ta exp[-|(x 1 ,x 2 ) r C'- 1 (xi,x 2 )] +0(e 2 ) 
(27r)^det(A) 



= ^ (xi) / u (x 2 )exp[exf A x aA x x 2 ] + 0(e 2 ) (A.5) 

where 

//.r-vl = 

2nVdetA 

is the marginal probability density of x = (0, 0') 



//o(x) = ^ — rexpf-^x 7 !"^] (A.6) 
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The correlation function (F(s)F(s')) = (^ r (xi)^ r (x 2 )) is 



(F(s)F(s')) = J d Xl j dx 2 ^( Xl ,x 2 )^( Xl )^(x 2 ) . (A7) 

The expectation (F(s)) is given by an analogous expression, in which ^(x) is integrated 
over x with weight /io( x )- Expanding the exponential in (A. 5), and using the fact that 
{F(s)) = gives 

(F(s)F(s')} =e J dx x J dx 2 ^ O (xi)^o(x 2 )^(x 1 )^(x 2 )xfl- I al- I x 2 + O(e 2 ) 

= eg T I- 1 cL4- 1 g + 0(e 2 ) (A7) 

where 

g= J dxxf(x)fi (x) . (A8) 

Equation (A. 7) shows that the correlation function {F(s)F(s')) is 0(e), implying that it 
is bounded by a multiple of the most slowly decaying correlation function of the fields 
occuring as arguments of the function T . It follows that the correlation function of F(s) 
will decay sufficiently rapidly at infinity if the correlation function of the fields has a 
sufficiently rapid decay. 

This approach extends directly to non-Gaussian fields in the commonly encountered 
case where the joint probability density factorises in the limit \s — s'\ — > oo, such that 
/i(x,x') =/i (x)//o(x')[l + 0(e)]. 
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